Escape time in anomalous diffusive media 
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We investigate the escape behavior of systems governed by the one-dimensional nonlinear diffusion 
equation 9ip — ^^^[fli^f/p] + 739^p'', where the potential of the drift, i7(a;), presents a double- well and 
D, V are real parameters. For systems close to the steady state we obtain an analytical expression of 
the mean first passage time, yielding a generalization of Arrhenius law. Analytical predictions are in 
very good agreement with numerical experiments performed through integration of the associated 
Ito-Langevin equation. For v ^ \ important anomalies are detected in comparison to the standard 
Brownian case. These results are compared to those obtained numerically for initial conditions far 
from the steady state. 
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I. INTRODUCTION 

The old problem of surmounting a potential barrier, 
known as Kramers' problem, is doubtlessly relevant in 
connection with many topics, in fields ranging from 
physics to finance. It is a key ingredient to understand- 
ing phase-transitions in complex systems, both in and 
far-from thermal equilibrium. In particular, the quantity 
known as the escape time (or mean first passage time) 
from one stable state to another has found numerous ap- 
plications in a variety of interesting and novel problems. 
For example, it plays a key role in stochastic resonance 
0, in describing fluctuation- induced transport such as 
occurs in kink motion |^ and ratchets Even the ex- 
tent of chaos in Hamiltonian systems has been shown to 
have connections with this quantity |^ . A nice collection 
of these and other stochastically driven processes can be 
found in Ref §. 

However, all of the above examples have been formu- 
lated within a standard Brownian framework, for which 
diffusion properties are normal. In this paper we look at 
the problem of calculating the escape time for systems 
exhibiting anomalous diffusion of the correlated type (in 
contrast to Levy type diffusion, which we do not dis- 
cuss here) . An understanding of escape time properties 
in such systems could open the door for understanding 
new stochastically driven phenomena. To our knowledge 
there has yet been little work done along these lines, al- 
though we are aware of some studies relating the anoma- 
lous transport properties on a random comb to the dis- 



tribution of mean first passage times Q . 

The systems we are interested in are such that the 
diffusion is dependent on the density of particles p, re- 
sulting in a diffusion coefficient which is proportional to 
a power (y — 1) of p. Many physical systems are well- 
described by this class of processes. Let us mention, 
amongst other examples, percolation of gases through 
porous media (z/ > 2) ||^, thin saturated regions in porous 
media {y — gravitational spreading of thin liquid 

films (y = 4) [||, heat transfer by Marshak waves [y = 7) 
|l0| , surface growth (z^ = 3) jOl, spatial diffusion of bi- 
ological populations {v > 2) [|12| , plasma flows {v < 1) 
|l3| . Explicitly, these processes are ruled by an equa- 
tion of the type known in the literature as porous media 
equation fl^ 



dtp{x,t) = Ddl[p{x,t)Y, 



(1) 



where a: is a dimensionless coordinate representing a 
bond-length, angle or any other chemical or physical state 
variable, t is the dimensionless time and vD > 0. Rewrit- 
ting the nonlinear term as dx{Dvp'^^^dxp), it becomes 
evident that the restriction Dv > guarantees that the 
ffux will be from more dense to less dense regions. 

Since the non-linearity in p is known to lead to anoma- 
lous diffusion if 7^ 1 (namely superdiffusion for < 1 
and subdiffusion for > 1 ||l5|,0, as < x^{t) >(x t^) 
important anomalies are also expected when crossing 
over a barrier is involved. Precisely, we want to unveil 
here how escape properties are altered when v ^ I. 
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The paper is organized as follows. In Sec. II we present 
the systems of interest and discuss some of their general 
features. Because fluctuations are determined by p{x, t) 
for ^ 1, the escape behavior will depend on the initial 
condition p{x,Q). Therefore we first consider systems in 
the vicinity of the steady state, a condition which allows 
analytical treatment. Numerical and analytical results 
for this case are presented in Sees. Ill and IV, respec- 
tively. In Sec. V we study numerically the escape be- 
havior of systems far from the steady state, comparing 
the results with the previous ones. Finally, section VI 
contains concluding remarks. 

II. THE SYSTEM 

Let us consider a set of identical particles immersed 
into a thermal environment such as that described by 
the porous media equation (|l|) . Under the influence of an 
external bistable potential U{x), introduced in order to 
probe the escape behavior, the density of particles evolves 
following the nonlinear Fokker-Planck (FP) equation: 

dtp{x,t) = d,[d^U{x)p{x,t)] -f Ddl[p{x,t)Y. (2) 

This class of equations has been the object of diverse 
previous studies [p^-|l7|. 

The stationary solution of Eq. (||) is 

Ps{x) = [l-{v-l)l3V{x)]t'' IZ, (3) 

where [/]+ — r7iaa;{/, 0}, Z is a (positive) normalization 
constant, P = Z''-^/{vD) and V{x) = U{x)-Uo, with Uo 
the absolute minimum of the potential. In the limit v — > 
1 the standard linear Fokker-Planck equation is obtained. 
In such a case, the steady state characterized by the 
Boltzmann-Gibbs distribution Ps{x) ^ exp(— C/(a;)/_D), 
is recovered. However, for v ^ I, the stationary solutions 
of Eq. (||) have the form of the Maximum Tsallis Entropy 
probability distributions, as already discussed previously 
|p5|-p7t, even in the absence of external drift [ p"4HT^ . It 
is worth recalling that phemonena such as full developed 
turbulence [ p^ , the hadronic transverse moment distri- 
bution in high energy scattering process e"'"e~ — > hadrons 
, among others, have been satisfactorily described in 
terms of distributions similar to (|^) instead of the canon- 
ical stationary one. 

Steady state solutions are illustrated in Fig. 1 for a 
quartic potential. Note that a cut-off condition (Tsallis 
cut-off), yielding regions with null probability, arises in 
the u > 1 case (see Fig. lb). For a quartic potential the 
condition v > —3 must hold so that the solutions can 
be normalized. However, the free-particle case requires 
J/ > — 1 so we restrict our discussion to this regime. 

The nonlinearity in the diffusion term of Eq. ^ ac- 
counts for the fact that the environment presents some 



kind of disorder or long range correlations in space- 
time leading to diffusion anomalies. The expression 
(3 — Z'^'^ / [vD) can be interpreted as a generalized 
Einstein relation for this scenario. Note that in disor- 
dered or correlated systems such as those discussed here, 
the standard Einstein relation is expected to be recov- 
ered in the absence of disorder pl[ |. This corresponds 
to the case of = 1 yielding the well-known result 
D = 1/(3. Also, as was shown in [Q, the time-dependent 
form of these Einstein relations can be used to demon- 
strate the anomalous scaling properties of these nonlin- 
ear diffusion systems. For the free particle one obtains 
< a;2(t) >(x l/Pit) cx Z'^{t) cx . 

The Ito-Langevin (IL) counterpart of Eq. (||) reads |l^] 

i = -a,C/(a;) + /^[p(x,t)]^77(i), (4) 

where r](t) is a delta correlated Gaussian noise with zero 
mean and variance 2. In the particular case i> = 1, the 
standard Langevin equation for constant noise is recov- 
ered. It is noteworthy that this is a phenomenological 
description, in which the microscopic trajectories are de- 
termined by the macroscopic quantity p when v ^ 1. 
Physically, this represents a kind of statistical feedback. 
As with state-dependent noise, it is to be seen as the 
influence of the environment, which is otherwise not ex- 
plicitly taken into account by the equations of motion. As 
a particle evolves, it interacts with the environment such 
that it reacts to the collective density of states around it. 
We can think of the subdiffusive case as a kind of " attrac- 
tion" to the other particles: Particles tend to stay close to 
the other particles, fluctuating not far from them. Con- 
versely, we can think of super-diffusive cases as a kind of 
reaction to the sparseness: If the particle is in a highly 
populated region then it is in a sense confined by the 
other particles, and fluctuations are not so large, but as 
soon as it gets into less dense regions it does not experi- 
ence this confinement and fluctuations can get very large. 

III. NUMERICAL RESULTS IN THE VICINITY 
OF THE STEADY STATE 

For numerical experiments we chose as prototype of 
double-well potential the quartic polynomial V{x) = 
ax'^ + bx^ + cx^ + d. The coefficients were chosen as 
in Fig. 1, for which {xL,xo,xii) = (0,1,3), with xl, xq 
and xr corresponding to the bottom of the left-hand well, 
the top of the barrier and the bottom of the right-hand 
well, respectively. We studied the escape behavior close 
to the steady state. That is, once a population of a large 
number of particles has already attained the steady state 
described by Eq. (H), a probe was injected a.t xl- Then 
its trajectory was obtained by solving, following the nu- 
merical scheme in ref. the IL (^ for p{x, t) = Ps{x), 
starting from x(t = 0) = xl- Typical trajectories are 
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displayed in Fig. 2. For v > \, fluctuations are reduced 
and trajectories result confined to the region within the 
cut-off boundaries (see also Fig. lb); moreover, when the 
diffusion constant D is smaller than a critical value D^. 
(here Dc — 0.17 for i/ = 2), the state space becomes dis- 
connected and crossings become forbidden. For v < 1, 
the amplitude of noise is enhanced in the regions of low 
density and the entire space tends to be populated. 

We measured the mean first passage time, i.e., the av- 
erage time interval T{xl — > x) that a particle at xi, 
takes to reach for the first time a given state x > x^. 
In Fig. 3 we present plots of T{x) = T{xl x) vs. x. 
For > 1 (Fig. 3a), plateaux become evident as D ap- 
proaches Dc indicating that most of the time is spent 
overcoming the barrier around xq- On the other hand, 
for V < \ (Fig. 3b), the passage time is sensitive to the 
exact final state and there is not a well defined plateau, 
even in the small-D regime. Moreover, as D decreases, 
the curves collapse to a limiting one for states below xr^ 
but grow faster above x/j, diverging in the limit I? — + 0. 
The escape behavior seems to be discontinuous at Z? = 0. 
In fact, for D = there is no diffusion, however, for finite 
D the particle is attracted towards the deepest valley at 
xn and becomes trapped within a typical time interval 
which is bounded from above. This effect can be un- 
derstood having in mind that fluctuations depend on D 
not only through the factor y^jl?! but also by means of 
the density through a factor that, iox v < 1, becomes 
very large outside the neighborhood of the absolute min- 
imum where particles tend to concentrate as I? — > 0. In 
other words, the deterministic case is not recovered when 
D ^ Q since the effective diffusion coefficient Dp'^~^ does 
not vanish in that limit due to the singularity at p = 0. 

IV. ANALYTICAL CONSIDERATIONS 

Let us show that these results can be understood ana- 
lytically. For a system in the vicinity of the steady state, 
we can consider the following approximation for Eq. (^ 

dtp{x,t) ~ d,[d^u{x)p{x,t)] + Ddl[{p,{x)Y-^p{x,t)]. 

(5) 

Once the FP equation is linear, the problem of escape 
from a well can be treated directly, following the same 
lines as for homogeneous processes characterized by time 
independent drift and diffusion coefficients . Basically 
an equation for the probability that the particle is still 
within a given interval of state space at time t is found 
using the corresponding backward Fokker-Planck equa- 
tion and solved under appropriate boundary conditions. 
In this way, one finds that the mean first passage time 
T{xi X2), for xi < X2, is given by 

T(xi^X2) = W\f3 f '[l-{iy-l)l3V{y)]f'dy 



X r [l-{u-l)pV{z)]f'dz, (6) 



where ^ = 1 if > and ^ = 1 — 2i^ if < 0. Expres- 
sion (^ reproduces numerical experiments with excellent 
agreement as illustrated in Fig. 3. 

In Fig. 4 we show T = T{xii) = T{xl Xr) as a 
function of 1/D (full lines), for different values of z/ > 0, 
as calculated from Eq. ^ . T represents a measure of the 
escape time from the left to the right-hand well, even in 
the V < 1 cases where plateaux are not well defined. In 
the range u > 1, T diverges at a value Dc^ defined by 
the cut-off prescription, below which the right-hand well 
becomes inaccessible. In the Q < v < \ case, T satu- 
rates as 1/D increases. The hyperdiffusive regime v <Q 
(hence D < Q), where spreading is faster than ballistic, 
demonstrates the same general features discussed for the 
region < < 1 but |D| must be considered instead of 
D. For any v and small 1/\D\ the escape time follows 
the power law T (3^ ^ 1/\D\'^ . 

li Xl ~ XL and X2 — xr, then, it is possible to find 
an approximate expression for the escape time T when 
1 1? I (hence 1//3) is sufficiently small, noting that the in- 
tegrands in Eq. (^) present sharp peaks at xq and xl 
respectively. In that case the integrals can be evaluated 
by a saddle-point approximation extending the integra- 
tion limits to the whole space. Following this procedure 
we arrive at 

^ui:JdW\ + ^^\i-{i^-i)f3v{xL)) 

where ul and luq are the frequencies at the bottom of 
the left well and at the top of the barrier, respectively. 
Expression (^ is a generalization of the Arrhenius law, 
which, as expected, is recovered in the limit i' 1. In 
fact, in that limit, T ~ {2n/ y^uJijIJo) exp{AV/ D), where 
AV = V{xo) — is the barrier height. 

For comparison, the approximation given by Eq. (|^) 
is also exhibited in Fig. 4 (dashed lines). The approx- 
imation is good for large l/l^?], as expected. It works 
better for v > 1. Let us comment the main features re- 
vealed by this expression. When v > 1, it foresees the 
divergence of T at finite D. In fact, Dc is obtained from 
l//3c ~ (i^ — l)V{xo)- When v < 1, saturation of T for 
large l/\D\ is also predicted (unless V{xl) =0) since /3 is 
an unbounded increasing function of 1/|D|. lfV{xL) = 0, 
then Eq. (^ indicates that T diverges for vanishing \D\. 

In particular, if < v < 1, T ^ /J^Ti^ ^ l/D^^ and 
the deterministic limit is achieved. In the limit v 1 the 
exponential growth of T with 1 /D is always recovered. 
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V. NUMERICAL RESULTS FAR FROM THE 
STEADY STATE 

The problem in the vicinity of the steady state ac- 
tually corresponds to a linear one with a state depen- 
dent diffusion coefficient. However, it allows an analyt- 
ical treatment which can be had in mind as a reference 
when studying more general cases. In order to test how 
the previous results compare to those of a more general 
situation, we also performed numerical studies of the es- 
cape properties far from the steady state. Particularly, 
we studied the case where particles are injected all at 
the same time at xl- This instance requires simultane- 
ous integration of the FP equation, in order to follow the 
evolution of t) starting from p{x, 0) = 5{x — xl), to- 
gether with integration of the IL equation (^), starting 
from x{t = 0) = XL = 0. Now, the parameter v must 
lie in the region > due to the divergence in Eq. (j^). 
An implicit finite-difference scheme with centered space 
differences was employed for numerical integration of the 
nonlinear FP equation [Q. The time evolution of the 
density is illustrated in Fig. 5. 

The escape time T as a function oi 1/D (symbols) ob- 
tained for different values of was included in Fig. 4. 
Let us compare this case to the precedent steady one. 
For sufficiently large D, T is not sensitively dependent 
on the initial distribution and Eq. (j6|) fits well to the 
numerical results for any v > 0, following the power law 
T ^ derived above. On the other hand, for small 

D, crossing times become closer to those of the standard 
case v ~ I for any ly. This can be understood as follows. 
For z/ > 1, passage times are smaller than those given 
by Eq. since, as the distribution evolves, there is an 
initial passage even between regions disconnected at the 
steady state (see Fig. 5(a)). However, our results suggest 
that the divergence of T for a finite critical D, close to 
Dc, still occurs. On the other hand, in the range i> < 1, 
crossing times are larger than those given by Eq. (|^) 
since now the density of particles is initially unfavorable 
for surmounting the barrier (see Fig. 5(b)). Saturation is 
not observed and the escape time increases with l/D ap- 
parently following a power-law. It is worth noting that, 
as derived above, a power-law with exponent 1/(1 — z^) is 
the one expected if the average effective potential felt by 
crossing particles has the absolute minimum at xl which 
is consistent with the observed density evolution (see Fig. 
5(b)). 



VI. FINAL REMARKS 

Summarizing, we have obtained the escape time for 
systems exhibiting anomalous diffusion due to a stochas- 
tic nonlinear dependency on the density. For steady-state 
conditions, we obtain an analytical expression for the 



mean first passage time whose predictions are in excellent 
agreement with numerical results (Fig. 3). This analyti- 
cal expression yields a generalization of Arrhenius law. A 
behavior quite different from that of the standard Brow- 
nian case v — 1 is depicted. Under close to stationary 
conditions, two regimes are detected: In the region v < 1 
(superdiffusion) , the escape time T saturates for vanish- 
ing D, if V{xl) 7^ 0, and grows with 1/D following a 
power-law, otherwise. In the region i> > 1 (subdiffusion) , 
T diverges at Dc (Fig. 4). 

These results give hints on what should be expected 
in more general cases. For systems far from the steady 
state, T grows with 1/D apparently following a power- 
law in the superdiffusive cases while T diverges at finite 
D in the subdiffusive ones. 
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CAPTIONS FOR FIGURES 

Figure 1: The cut-off condition, (a) Dimensionless 
double- well potential V{x) = ax* + bx^ + cx^ + d, with 
a = 1/48,6 = -1/9, c = 1/8, d 3/16. The stationary 
distribution Ps{x) is shown for = 2 (b) and 0.5 (c), 
for different values of D as indicated on the figure. For 

< 1 the full state space is covered with power-law tails. 
For ly > 1 a cut-off restricts the attainable space. Ob- 
serve in (b) that as D decreases particles become more 
confined until only the neighborhood of the deepest val- 
ley is allowed. The horizontal lines in (a) represent the 
cut-off condition V{x) = 1/(3 which defines the allowed 
regions for i/ = 2 and the same values of D as in (b) . All 
quantities are dimensionless. 

Figure 2: Typical trajectories x vs. t for {ly, D) = 
(0.5,0.5) (dark gray), (2,0.5) (black) and (2,0.15) (light 
gray). 

Figure 3: T{x) = T{xl x) vs. x for different val- 
ues of D indicated on the figure and = 2 (a), 0.5 (b). 
Circles correspond to numerical experiments (mean value 
over 1000 realizations) and full lines to theoretical pre- 
diction given by Eq. (^). 

Figure 4: Escape time T = T{xii) sls a function of 
1/D, for different values of z^ > indicated on the fig- 
ure. Full lines are generated from Eq. (|^). Dashed lines 
correspond to the low- 1? approximation given by Eq. (j^). 
Symbols correspond to the initial condition where all the 
particles (at least 1000) are injected at the same time at 
Xl- Dotted lines are guides for symbols. Insert: Detail 
(semi- log) of the \ow-D region for v <1. 

Figure 5: Time evolution of the density of parti- 
cles obtained by numerical integration of Eq. (^ with 
p{x,0) = 5{x) for (i/,D)=(4.0,2.5) (a) and (0.5,0.1) (b). 
The profiles correspond to times t indicated on the figure. 
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